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An analysis is given of the Stueckelberg extension of the Standard Model with a hidden sector 
gauge group C/(l)x where the mass growth for the extra gauge boson occurs via the Stueckelberg 
mechanism, and where the kinetic mixing in the U{l)x x U{1)y sector is included. Such a kinetic 
mixing is generic in a broad class of supergravity and string models. We carry out a detailed global 
fit of the model with the precision LEP data on and off the Z pole, with within 1% of the of 
the Standard Model fit. Further, it is shown that in the absence of matter in the hidden sector, there 
is a single effective parameter that controls the deviations from the Standard Model predictions, and 
the dependence on the kinetic mixing emerges only when matter in the hidden sector is included. An 
analysis is also given of milli-charged dark matter arising from the hidden sector, where it is shown 
that such dark matter from the Stueckelberg extension can satisfy WMAP-3 data while allowing for 
a sharp Z' resonance which can be detected at the Tevatron and at the LffC via a dilepton signal 
generated by the Drell-Yan process. 



I. INTRODUCTION 

Recent works on the Stueckelbergextension of the 
SM [H (StSM) and of the MSSM 0, i] (StMSSM), 
have shown consistency with the LEP data while al- 
lowing for the possibility of a narrow resonance which 
could lie as low as just above the Z mass Jj, iJ, 0, [H . 
Similar phenomena regarding a narrow resonance are 
seen in other classes of models such as those based on 
universal extra dimensions [6], models with a shadow 
sector 0, H, Q , and in the models considered in Ref . 
[13 . The possibility of a narrow graviton resonance 
also arises in the RS model [ll|, [l^ • Thus the study 
of narrow resonances is a topic of significant inter- 
est. In this work we focus on the Stueckelberg ex- 
tensions further. Such models may be a field the- 
ory realization of string, models arising from orien- 
tifolds [H, Hi, [H, [H, |l3 and their discovery could 
be harbingers of a new regime of physics altogether. 
The supersymmetric Stueckelberg extensions have the 
possibility of generating a new type of neutral Majo- 
rana fermion which is massive and extra-weakly inter- 
acting, and with R parity conservation, a candidate 
for cold dark matter [l3|. Indeed a detailed analysis 
shows that the relic density in such models is con- 
sistent with the three year WMAP data Addi- 
tionally, the Stueckelberg extensions give rise to milli- 
charges for matter residing in the hidden sector [JQ. 
A recent analysis [20,] indicates that such matter can 
annihilate in sufficient amounts to satisfy the relic 
density constraints from WMAP-3 for a broad Z' res- 
onance. 

The main focus of this work is an extension of the 
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class of models considered in Refs. [l|, d, H, 0, Q by 
including kinetic mixing between the two Abelian 
U{1) gauge fields. Specifically we consider the 
extended electroweak sector with the gauge groups 
SU{2)l X U{1)y X U{l)x where the Stueckelberg 
mechanism along with the spontaneous breaking in 
the Higgs sector generates the vector boson mass, 
and a mixing in the gauge kinetic energy of the 
U{l)x X U{1)y sector is included. Such kinetic 
mixings can arise in a variety of ways pll . [23 | and 
enter prominently in models of the type considered in 
Refs. 8, 9]. The model considered here encompasses 
the models of Refs.P, i, i, H which can be 
obtained in certain limits of the model discussed here. 
Inclusion of the kinetic mixing in the Stueckelberg 
extension enhances significantly the parameter space 
where new physics can exist consistent with the 
stringent LEP, Tevatron, and WMAP constraints. 
This parameter space includes the possibility of a 
narrow Z' resonance ver y d istinct from the Z' of the 
conventional models [H,]!!!!!. 

The outline of the rest of the paper is as follows: In 
Sec. (II) we give a description of the Stueckelberg ex- 
tension with kinetic mixing in the U{l)x x U{1)y sec- 
tor. In Sec. (Ill) we give a detailed numerical analysis 
of the electroweak constraints from LEPI and LEPII. 
In Sec. (IV) an analysis of milli-charged dark matter 
that arises from the hidden sector is given. Conclu- 
sions are given in Sec.(V). Appendices contain several 
mathematical details, including an analysis regarding 
the origin of milli-charged matter. 

II. THE STUECKELBERG Z' EXTENSION 
WITH KINETIC MIXING 

In this section we discuss the U{l)x Stueckelberg 
[26j extension of the Standard Model (SM) with gauge 
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kinetic mixing (StkSM). We assume that the quarks, 
leptons and the Higgs field of the SM do not carry 
U{l)x quantum numbers, and the fields in the hidden 
sector do not carry quantum numbers of the SM gauge 
group. Thus the U{l)x sector is hidden except for 
the U{l)x X U{1)y mixings in the gauge vector boson 
sector as given by the effective Lagrangian 
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Here i?^ is the U{1)y gauge field, is the U{l)x 
gauge field with coupling strength gx to hidden mat- 
ter through the source Jx, o is the axion, M\ and 
M2 are mass parameters, and 8 is the kinetic mixing 
parameter. The a field is charged under both U(X)x 
and f/(l)y and the Lagrangian of Eq. (1) is invariant 
under the U{V)x x C/(l)y gauge transformations 

8B^, = df,\x, SCf, = 0, S<j = -AhXx 

6Bf, ^ 0, SC^ = df^Xy, 5a ^ -MiAy. (2) 



which diagonalizes the kinetic terms has the form 
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1 



c, = 



1 



Ss = SCs. (5) 



The matrix R is then defined by the diagonalization 
of the mass matrix 
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The model of Eq.(l) involves three parameters: 
Mi,M2,<5 or alternately Mi,e, (5 where e = M2/M1. 
Expressing the transformation T KR in terms of 
the matrix elements, one has 



Tij = CsRij, 

T2j ~ R2j — SsRij, 



(7) 



The neutral current interaction with the visible sector 
fermions is given by 



Cnc = --Y^if^Df^rh + {L^ R)], 



(8) 



The above model has a non-diagonal kinetic mixing 
matrix (/C) and a non-diagonal mass matrix (Mg^) and 
in the unitary gauge [l| in the basis — (C, B,A^) 
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A simultaneous diagonalization of the kinetic energy 
and of the mass matrix can be obtained by a transfor- 
mation T = KR, which is a combination of a GL{3) 
transformation (K) and an orthogonal transformation 
(R). This allows one to work in the diagonal ba- 
sis, denoted by E where E'^ = {Z',Z,A^), through 
the transformation V = {KR)E, where the matrix K 



where is the covariant derivative with respect to 
SU{2)l X U{1)y X U{l)x gauge group except that, as 
mentioned in the beginning, we assume that the visi- 
ble sector matter, i.e., quarks, leptons and the Higgs, 
are not charged under U{l)x- Thus the covariant 
derivative includes only the SU{2)l gauge coupling 
g2 and U{1)y gauge coupling gy- The diagonalization 
also leads to the following relation for the electronic 
charge 
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Thus gy is related to by 
gy = 7V1 + < 



2Se, 



(9) 



(10) 



and one may write the neutral current interaction so 
that 



J 



Cnc = ^^^^Ir [Wf - l5a'f)Z'^ + {vf - j5af)Z^] f + ef^QfA^J, 



(11) 
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gyT22)Tf + 2gyT22Qf], 

9YT22)Tf], 

9YT2i)Tf + 2gyT2iQf], 



:[i92T31 - gyT2i)Tf], 



(12) 



and where, as usual, — Tj + Yf/2. We note that 
the mass and kinetic mixing parameters enter not only 
through Tij but also through gy via the constraint of 
Eq.®. 
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Z Mass,1% LEPI and LEPII Constraints [ m1 = 200 GeV ] 




FIG. 1: An analysis of in StkSM model in the e - 5 
plane. The center white region is where of the StkSM 
model is within 1% of the SM fits. Along the line e = 5 
the Z' decouples and the model gives the same fit to 
data as in the SM (see also Appendix A). 
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FIG. 2: The plots give an analysis of the LEPII constraint. 
The upper plot, which has a Z' mass range of .2 — 1 TeV, 
shows the relevant LEPII contact interaction parameter 
A^y as a function of for the 19 observable of Table(l), 
where a ^ 25 is the SM fit as given in Table(l), and 
where the (yellow,red) shaded regions correspond to A^'^ 
= (21.7, 17.1) TeV [131. The lower plot for in the range 
(25-50) gives A^'^ as a function of the Z prime mass. 



III. CONSTRAINTS FROM 
ELECTROWEAK DATA 

We discuss now the constraints on the model 
with both mass mixing and kinetic mixing from the 
precision electroweak data. We start by assuming 
that the hidden sector does not contain matter, and 
the case when matter is included in the hidden sector 
is discussed in Sec. (IV) and in Appendix B. To obtain 
the allowed range of e and <5, wc follow the same 
approach as in Ref. 0, [23|. The first constraint 
comes from the comparison of the one sigma error in 
the prediction of the Z boson mass in the Standard 
Model and a comparison of this result with experi- 
ment, leads to an error corridor, 5Mz ^ 37 MeV, 
where one can accommodate new physics. However, 
the more stringent constraint comes from fits to the 
high precision LEP data on the branching ratios of 
the Z decay and from the various asymmetries at 
the Z pole, when one demands that the x^ fits of 
StkSM are within 1% of that of the Standard Model. 
We will refer to this as the LEPI 1% constraint 
in the rest of the paper. Details of the method 
employed for the electroweak fits can be found in 
the analysis of Ref. [1, and here we present just 
the results in the extended model with kinetic mixing. 

Table 1 gives a fit to the LEP data for a spe- 
cific point in the Stueckelberg parameter space with 



e = .06, (5 = .03, and Mz' « Mi = 200 GeV. In the 
analysis we have taken into account the constraint be- 
tween gy and gy^^ and the inclusion of this constraint 
improves the electroweak fits over that of previous 
analyses for the case 5 = ^, "5^. Thus the analy- 
sis of Table 1 shows that in the StkSM one finds x^ 
fits which are at the same level as in the SM. An anal- 
ysis of x^ in the LEPI fits in the e — 5 parameter space 
is given in Fig.([T]). Specifically Fig.([T|) shows that a 
large region of the parameter space can satisfy the 
LEPI 1% constraint. A striking aspect of the analy- 
sis of Fig.(IT|) is that this constraint is satisfied even 
though e and 6 can get significantly large, as long as 
(e — (5) is small. The physics of this is explained in 
Appendix A where it is shown that in the absence of 
matter in the hidden sector, there is only one effective 
parameter, e, that enters the analysis of electroweak 
physics. 

We discuss next the LEPII contraints. These con- 
sraints are typically characterized by the parameter of 
contact interaction A, and the LEPII group finds that 
t^vv > (21.7, 17.1) TeV ^ to be the most constrain- 
ing. The StkSM model predicts the theoretical value 
of Avv through the following formula 
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TABLE L Fits to 19 Z pole observables. Column 2 is given by the PDG [28|, while the data in column 3 is from the SM 
Fit of the LEP EWWG ^S]- The column labeled St Fit is an analysis for the input e = 0.06, & = 0.03, and i\fi = 200 
GeV. In the last column PULL is defined by (Experiment - FIT)/A, and = E PULL^. 
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A numerical analysis of the LEPII constraints is given 
in Fig. HI) and Fig.®. The analysis of Fig. ([1]) exhibits 
that the LEPI 1% constraint is more stringent than 
the LEPII constraint, and thus the LEPII constraint 
is automatically satisfied once the LEPI 1% constraint 
is satisfied. This result is supported by the analysis of 
Fig-© which shows that the value of Ayy predicted 
by the model in the parameter space consistent with 
the LEPI 1% constraint is significantly larger than 
the lower limit of the LEPII constraint. The blue 
points that enter the shaded regions are eliminated 
by the LEPII constraint. However, these points also 
correspond to large fits to the LEPI analysis and 
are eliminated by LEPI 1% constraint as well. Thus, 
for a narrow Z' , the LEPI 1% constraint is stronger 
than the LEPII constraint. 



IV. MILLI-CHARGED DARK MATTER 
FROM THE HIDDEN SECTOR 

In the previous section we did not include matter in 
the hidden sector which is defined as matter which is 
neutral under the SM gauge group but carries U{l)x 
quantum numbers and thus couples only to C^. The 
kinetic and mass mixings in the U{l)x x U sectors 
typically generate milli-charges for such matter. The 
conditions for the origin of milli-charges arising from 



such mixings are discussed in Appendix B, where sim- 
ple examples are worked out to explain the constraints 
that lead to the appearance of such charges. Milli- 
charges have been examined in many works both the- 
oretically and ex p erimentall y [t I -IsiIIs^ . [H, HI, [H, [H, 
m m, m SO, mi m, S, Hlllla Most of these 
analyses are in the context of kinetic mixing model of 
\J\. Here we consider the milli-charged matter in the 
hidden sector within the context of the Stueckelberg 
extension of the SM with both mass and kinetic mix- 
ing. If milli-charged matter exists then both the Z and 
the Z' can decay into it if kinematically allowed to do 
so. For the mass scales we investigate the milli-charge 
particle has a mass larger than M2/2. In this case all 
of the electroweak constraints discussed in Sec. HI are 
unaffected. Further, the Z prime can decay into the 
millicharged matter if the mass of the hidden matter 
is less than Mz' /2. Such decays increase the Z' width 
and thus decrease the branching ratios of the Z' decay 
into the visible sector which depletes the dilepton sig- 
nal in the Drell-Yan process. A relatively strong dilep- 
ton signal manifests in the analysis of Refs. [3, 0, 13| 
where the Z' decays into the hidden sector were taken 
to be comparable to the Z' decays into the visible sec- 
tor, i.e., ~ rp?. This constraint then leads to a 
sharp Z' resonance, but the decay of the hidden sector 
matter via the Z' pole is not strong enough to anni- 
hilate the hidden sector milli-charged dark matter in 
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FIG. 3: An analysis of the relic density of milli-charged 
particles arising in the StSM Z' model from the hidden sec- 
tor for the case 5 = 0, M;^ = 150 GeV, e = (.01 - .06) with 
.01 for the innermost curve and moving outward in steps of 
.01. The (yellow, grey) regions (M^/ < 2M^, Mz' > ^M^) 
correspond to a (narrow, broad) Z' resonance, and the 
WMAP-3 rehc density constraints are satisfied for both 
a broad Z' resonance and a narrow Z' resonance as ex- 
hibited by the 2a WMAP-3 red and black bands. The 
region of narrow Z' resonance is constrained by the LEP 
and Tevatron data. The region in the 2a WMAP-3 band 
can be probed via a dilepton signal as shown in Fig. ([4]). 
The red band to the left is excluded by the CDF 95% C.L. 
[4^ data while the black band is consistent or on the edge 
thereof, with all constraints which can produce an observ- 
able dilepton signal. 



sufficient amounts to be consistent with experiment, 
unless extreme fine tuning is used (we return to this 
issue at the end of this section). 

The recent work of Ref. [l^ has carried out an 
explicit analysis of putting a pair of Dirac fermions in 
the hidden sector, and made the interesting observa- 
tion that for values gxQx < 0(1) the decay width of 
Z' into the hidden sector Dirac fermions (x) can be 
of GcV size, and consequently the hidden matter can 
annihilate in sufficient amounts to satisfy the relic 
density. We have carried out a similar analysis using 
the thermal averaging procedure in the computation 
of the relic density as described in Appendix C. Our 
conclusions are in agreement with the analysis of Ref. 
[201] in the region of the parameter space investigated 
in Ref. j20| when no kinetic mixing is assumed in the 
absence of thermal averaging. In our work we take 
the kinetic mixing into account in the analysis of the 
relic density. We also make a further observation 
that there exists a significant region of the parameter 
space where it is possible to satisfy the relic density 
constraints and still have a narrow Z' resonance 
which can be detected at the Tevatron and at the 
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FIG. 4: An exhibition of the dilepton signal a ■ Br{Z' — » 
at the Tevatron consistent with the WMAP-3 relic 
density constraint as a function of Mz' when 2M^ = 300 
GeV. The curves in ascending order are for values of e in 
the range (0.01 — 0.06) in steps of 0.01. The dilepton signal 
has a dramatic fall as Mz' crosses the point 2M^ = 300 
GeV where the Z' decay into the hidden sector fermions is 
kinematically allowed, widening enormously the Z' decay 
width. The green shaded regions are where the WMAP- 
3 relic density constraints are satisfied for the case when 
there is no kinetic mixing. Red and blue regions are for 
the case when kinetic mixing is included. The D0 data set 
[49II was collected in the search for narrow resonances (RS) 
and is a stronger constraint to apply on this model than 
the recent CDF data which put constraints on the 
parameter space when the Z' can decay only into matter 
in the visible sector. 



LHC using the dilepton signal via a Drell-Yan process. 

We give now further details of our relic density anal- 
ysis. In the numerical analysis we will use gx = 9y^^ 
and Qx — 1 unless stated otherwise. We begin by 
considering the case when 5 = which is the StSM 
case. Fig. ([3]) gives an analysis of the relic density 
as a function of Mz' for the case = 150 GeV, 
and e in the range (0.01 — 0.06). Here one finds that 
the relic density is satisfied on two branches, one for 
Mz' > 2M^, and the other for Mz' < 2M^. In the 
region Mz' > the relic density is satisfied over 

a broad range of Z' masses for appropriate e values. 
In this region the decay width of the Z' is large, and 
thus the branching ratio into the visible sector is sup- 
pressed, which would make the dilepton signal from 
the Z' decay difficult to observe. This is exhibited in 
Fig. (HI) where one finds that the dilepton signal essen- 
tially disappears in the region Mz' > 2M^. Next we 
examine the region Mz' < 2M^. Here one finds once 
again a satisfaction of the relic density even though 
the Z' mass lies below 2M^. The very sizable range 
of the Z' mass over which the relic density can be 
satisfied arises due to thermal averaging over the Z' 
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pole. As expected in this region the Z' will appear 
as a sharp narrow resonance since the decay of the Z' 
into the hidden sector fermions is kinematically dis- 
allowed. An analysis of the dilepton signal for this 
case is also given in Fig. ([4]), and a more detailed view 
is given in Fig. ([5]), which shows that the Drell-Yan 
signal pp Z' e~^e~ is enormously enhanced for 
Mz' < 2M^. Thus we have a region here of the pa- 
rameter space where one will have a sharp resonance 
giving a visible dilepton signal while at the same time 
producing milli-chargc dark matter consistent with 
WMAP-3. 

We carry out a similar analysis for the case when the 
kinetic mixing parameter 6 is non- vanishing. Fig.([n]) 
gives an analysis of the effect of 5 on the thermally 
averaged cross section J(xf), defined in Appendix C, 
as a function of the Z' mass. One finds that J{xf) is 
affected in a significant way by the change in 5 and the 
particulars of the modification due to 5 are explained 
in Appendix D. In Fig. ([7]) an analysis of relic density 
as a function of AIz' is given for e = 0.04 and 5 in 
the range 5 — (.05 — .25). The analysis shows the 
relic density exhibits a strong 5 dependence. Quite 
significantly, the region in Mz' over which the relic 
density can fall in the WMAP-3 region is widened. 
Once again, in the region Mz' > 2M^ the dilepton 
signal will be too dilute to be observable as shown 
in Fig. ([7]). On the other hand, in the region Mz' < 
2M^, the Z' is a sharp narrow resonance since the Z' 
decay into the hidden sector fermions is kinematically 
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FIG. 5: Same as Fig. ([4} except that only the mass region 
with a detectable dilepton signal at the Tevatron is exhib- 
ited but additional S values are included in the analysis. 
The plots show that the region allowed by WMAP-3 con- 
straints moves to the right for positive 5 and to the left 
for negative S. This phenomenon is explained in Appendix 
D. The inclusion of kinetic mixing is seen to enlarge the 
parameter space where the relic density constraints are 
satisfied and where an observable dilepton signal at the 
Tevatron can occur. 



disallowed. Thus the dilepton signal here is strong 
as can be seen in Fig.(|l]) and Fig.® and should be 
observable with sufficient luminosity. 

An analysis for the case of the LHC in the Drell- 
Yan process pp ^ Z' ^ e+e~, is given in Fig.®. 
Here it is shown that the dilepton signal in the region 
Mz' < consistent with the WMAP-3 constraints 
would be discernible with 10/6^^ of integrated lumi- 
nosity for e > .04. The dilepton signal has a dramatic 
fall as Mz' crosses the point AIz' = 2M^ beyond 
which the Z' decay into the hidden sector fermions 
is kinematically allowed, widening enormously the Z' 
decay width. The green shaded regions are where 
the WMAP-3 relic density constraints are satisfied 
for the case when there is no kinetic mixing. Red 
and blue regions are for the case when kinetic mixing 
is included. When computing the dilepton signal, 
a 20 GeV mass bin is used, and a 50% detector 
cut is assumed. The thick solid line corresponds to 
the LHC discovery limit with 10/6^^ of integrated 
luminosity expected in the first year run, and the 
criteria used for the discovery limit is 5V A^s m or 10 
events, whichever is larger, where Nsm is the number 
of SM background events. The analysis shows that 
data in the first year running at the LHC should 
be able to constrain the parameter space down to 
e - 0.04 for = 500 GeV when Mz' < 1 TeV. 

An interesting issue concerns the question regard- 
ing how small e can be for WMAP-3 relic density 
constraints to be satisfied. Fig.® addresses this 
question where an analysis of the relic density is given 
when e = 10~^. One finds satisfaction of the relic 
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FIG. 6; Illustration of the variation with S of the ther- 
mally averaged cross section for the annihilation of milli- 
charged particles when Mx = 250 GeV and e = .03. The 
generic effect of 5 is to modify the shape of J{xf) and the 
region where the relic density constraints can be satisfied 
is enlarged. Further details are given in appendix D. 
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FIG. 7: An analysis of the relic density of milli-charged 
particles for the case when kinetic mixing is included in the 
Stueckelberg Z' model. The analysis is done for — 150 
GeV, e = .04, and d = (.05, .075, .10, .15, .20, .25), where 
the values are in descending order for AIz' > 300 GeV. The 
red and black bands are the WMAP-3 constraints where 
the black band also produces an observable dilepton signal. 
The analysis shows that for e fixed, increasing 5 increases 
the parameter space where the WMAP-3 relic density con- 
straint is satisfied, while allowing for a detectable Z prime 
signal as shown in Fie:s. (|4l5|) . 

density in this case, and even smaller e were found 
admissible. Further, while the dilepton signal at the 
Tevatron in this case will be suppressed, it could still 
be visible at the LHC with sufficient luminosity. 



Finally, we note that within the context of the 
Stueckelberg model it is possible to place indirect 
limits on the milli-charge coupling of hidden sector 
fermion with the photon from the Tevatron data. An 
analysis of the limits on the Stueckelberg mixing pa- 
rameter e was presented in [3] for the case ^ = 0. 
For this case, the milli-charge QmiiH, where Qmiine 
is the coupling of the photon with the hidden sec- 
tor fermions (see Appendix D) is determined to be: 
Qmiiii ~ £■ Thus one may directly translate the limits 
obtained in Q to limits on the milli-charged coupling 
of the hidden fermion with the photon. In the con- 
text of the present analysis, the cross-section predic- 
tions given here, (as for example, in Fig. ([5])) with their 
overlapping WMAP-3 bands, shows this explicitly. 



V. CONCLUSION 

In the above we have given an analysis of the 
Stueckelberg extension of the Standard Model with 



10^ 




900 950 1000 1050 1100 1150 



Z' Mass [GeV] 

FIG. 8: An exhibition of the dilepton signal as given by 
the number of events for 10/&^^ of integrated luminosity 
at the LHC consistent with the WMAP-3 relic density 
constraint as a function of Mz' when = 500 GeV. The 
curves in ascending order are for values of e in the range 
(0.01 — 0.06) in steps of 0.01. Criteria set for the discovery 
limit are discussed in the text. 
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FIG. 9: An analysis of the relic density of milli-charged 
particles for the case of a very small e. One finds that 
even for this case one can satisfy the WMAP-3 relic density 
constraints. In this case the hidden sector fermion mass 
is predicted to be essentially Mz' /2. While the dilepton 
signal from the Drell-Yan process is too weak to be observ- 
able at the Tevatron, it is strong enough to be detectable 
at the LHC with sufficient integrated luminosity 



inclusion of the kinetic energy mixing in the U{l)x x 
U{1)y sectors. Such kinetic mixings are quite generic 
in models with more than one L''(l) gauge group. It 
is shown that in the model with both the mass and 
kinetic mixing and in the absence of matter in the 
hidden sector, the sensitive parameter which measures 
the deviation from the Standard Model is given by e 
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as defined by Eq. ([22]) which is a specific combina- 
tion of e and d, where e measures the mass mixing 
and S measures the kinetic mixing. However, when 
matter in the hidden sector is taken into account, elec- 
troweak physics depends on both e and S. An analysis 
of the relic density of milli-charged dark matter which 
is generic in Stueckelberg extensions is given. Here our 
analysis is in agreement with the work of Rcf. 20] for 
the case when no kinetic mixing is taken into account. 
Inclusion of the kinetic mixing is seen to increase the 
region of the parameter space where the relic density 
constraints consistent with WMAP-3 can be satisfied. 
We also analyze the Z' signal. As noted in Ref. [2^ on 
the branch Mz' > the dilepton signal from the 
Z' decay is too small to be observed at colliders, and 
our results are in agreement with this analysis. How- 
ever, we note that on the branch Mz' < 2Af^, there 
is a significant region of the parameter space where 



the relic density constraints can be satisfied and the 
dilepton signal from the Z' decay via the Drell-Yan 
process is strong enough to be observed at the Teva- 
tron and at the LHC. The analysis also shows that 
relic density constraints can be satisfied for values of 
e as low as 10^^ and even smaller values are possible. 
An interesting issue concerns the detectability of such 
dark matter in laboratory experiments which might 
put further limits on the parameter space or on the 
component of the relic density such matter can con- 
stitute. However, such an investigation is outside the 
scope of this work. 
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I 

Appendix A: Details of the Stueckelberg extension with kinetic mixing 



In this Appendix we give further details of the Stueckelberg extension with kinetic mixing. The mass matrix 
in the vector boson sector, after applying the transformation that diagonalizes the kinetic energy, is given by 



= K'^Ml;tK, (14) 

and we display its explicit form below 



5) 



4 ' ^ 



1 2 5v2E2iHiI 



l92lV 



2 <5\/l-2£(5+e2 



(15) 



2e5 



492^ 



The eigenvalues of Eq. (fT5|) are 

M2 = 0, Af| ^{q- p)/2. Ml = {q + p)/2, (16) 
where 



P 



Ml [3 - 



^M2U_92^2p 

4. 



1^ 13 + 91 2 



(17) 



and where (3 = [l - 2t6 + e^) / {I - 6'^) . As expected in 
the limit e = = (5 the Stueckelberg sector decouples 
from the SM sector and one gets 



Mz = ^^^^ ^ V, Mz' = Ml. 



(18) 



The expression for Mz of Eq. (|T8l) is exactly the re- 
sult in the SM. However, quite remarkably decoupling 
takes place also in the limit when e and S are non- 
vanishing but satisfy the relation e = 6. Here, while 



the e, S dependence persists in the A/^ mass matrix of 
Eq. (fTS)) . the eigenvalues are still given by Eq. ((T8)) 
in this limit. Further, the coupling of the Z boson to 
quarks and leptons reduces to that of the SM in this 
limit. One may understand this result by examining 
the explicit form of the transformation matrix of Eq. 
(O which for the case e = 6 becomes 



/I 



92 



T : 



, T 1 

— 



92 1 

7 



\ 



J 



(19) 



A direct substitution in the Z and Z' couplings 
then shows that they are identical to the case of the 
SM. These analytic forms explain the limits seen 
numerically in the analysis of Fig.((T]). 

Next we consider the case with arbitrary mass mix- 
ing and kinetic mixing. Here before diagonalizing the 
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mass matrix it is convenient to perform the following 
orthogonal transformation (O) given by 







O 



6 VT~5^ 
1 



Mil 



(20) 



which transforms the mass matrix to = O'^ APO 



( Mf 

Mil Mil' + 2^7^(1 + e') -^527^1+1^ 



V 



-92 



4 92 



(21) 



where e is defined so that 



(22) 



We note that the mass matrix Eq. ((2T|) looks exactly 
the same as for the mass matrix one has if there was 
just the Stueckelberg mass mixing except that e is 
replaced by I. Thus the orthogonal matrix which di- 
agonalizes Eq. (|2ip will be the same as for the StSM 
case except that the free parameter is changed from e 
to e. We parameterize the rotation matrix TZ defined 
by n^M'^n = Diag{M^,,M^, 0) by the following 



n 



cos ip COS (j) — sin 9 sin (j) sin ^ sin ^ cos (f) + sin 6 sin (j) cos tp — cos 9 sin < 
cos ■0 sin + sin 9 cos sin ■0 sin ■0 sin — sin 9 cos cos V' cos cos (j) 
— cos sin tp cos cos -0 sin 9 



where the angles are defined so that 



tan 9 = 



92 



ta,n(j) — e, tan 2-0 = 



2sin6lMo2e 



Ml - M^ + {Ml 



Mi 



M,i 



(23) 



(24) 



and Mq = Mz{i ~ 5) ~ g\ + 7^/2, My/ = 92v/2. The transformation relating the initial basis and the final 
diagonal basis is ^ = [K{5)0{5)'R{l)]E, where = {C,B,A^), and E"^ = {Z',Z,A^). The neutral current 
interaction can now be written in the form 



CNc^J^sii,s)nii)E 

where = (gx-Tx, 7>/y , ff2>^2 )i ^^'^ given by 



S{l,S) = 







(25) 



(26) 



When Jx — 0, one finds that the neutral current interaction of Eq. l(25|l has no dependence on S. Expressing 
the tree level interaction in terms of the reduced vector and axial vector couplings as defined in Eq. (jlip. we 
find 



a/ 
v' 



cos-0 [(l-esin6'tanV;)r^-2sin^0(l-ecsc6'tanV')(5/] , 
cos -0 [1 — e sin 9 tan -ip] T^, 

-cos^p [(tanV- + esin6')T^ - 2 sin^ 6* (ecsc 6* + tan0) Q/] , 



a'f = — cos ip[ta,mp + € sin 9] T? . 



(27) 
(28) 
(29) 
(30) 



Thus in the absence of the hidden sector matter, the effect on the neutral current sector of including the kinetic 
mixing in the Stueckelberg model, is the replacement of e by e, and in the limit e ^ the Z couplings of the 
SM are recovered. However, with the inclusion of matter in the hidden, i.e., Jx 7^ 0, the neutral current sector 
will depend on both I and 6. We exhibit this below 



Cnc = C-Nc{Jx = 0) + 9xJx 



(7^ll - 5^7^2l) + (7^l2 - 5'^7^22) z^ + (7^l3 - 557^23) ai 



(31) 



where Ss appears in Eg. pip , and is as defined in on both e and 5 and an analysis of the branching 
Sec. (II). In this case the decay modes of Z' depend 
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ratios should reveal the presence of kinetic mixing. 

Further, as discussed in Sec. (IV) and in Appendix 
D, the relic density analysis depends sensitively on S. 



and leptons, while A' is the orthogonal massless 
state, and J' is the source in the hidden sector. Here 
the coupling of the photon with the hidden sector 
is proportional to d and thus the hidden sector is 
milli-charged if S is small. 



Appendix B: On the origin of 
milli-charged matter 

In this Appendix we illustrate the mechanism which 
generates the milli-charge in the context of the anal- 
ysis of this paper. We start with the kinetic mixing 
model 01 with two gauge fields Ai^,^2ai correspond- 
ing to the gauge groups U{1) and 11(1)'. We choose 
the following Lagrangian C ~ Co + Ci where 



Co 



- Fi F^" 



(32) 



To put the kinetic energy term in its canonical form, 
one may use the transformation 



At^ 



Kn 



A'f" 



(33) 



However, the transformation that canonically diago- 
nalizes the kinetic energy is not unique. Thus, for 
example, K = KqR instead of Kq would do as well 
where R is an orthogonal matrix 



R = 



Here £1 is given by 



cos W — sm ( 
sin 9 cos 6 



(34) 



£1 = A'^' 
+A^' 



cos 9 , 



',95 



smW — 



sin 95 



J' + cos 6* + I J 



.(35) 



In this case we see that each of the massless states 
interacts with the sources J and J'. However, one may 
choose 9 to get asymmetric solutions. For instance for 
the case 9 = arctan — <^^J one has 



£1 = 



+ A'''J' (36) 



In this case while A' interacts only with the source 
J', A interacts with both J and J', with the coupling 
to the source J' proportional to the kinetic mixing 
parameter 5. We identify A with the physical photon 
field, J with the physical source arising from quarks 



Next we consider a model with kinetic mixing where 
a Stueckelberg mechanism generates a mass term of 
the type considered in Eq. (1) 



C 



Mass = -hd^A^^A^; - ^M^A2^.A^ - MiAhAi^At^. 

(37) 

In this case diagonalizaton of the mass matrix fixes 9 
so that 



9 = arctan 



l-5e 



(38) 



and the interaction Lagrangian is given by 



1 



2^e + e2 VVT~^ 
1 



Ju 



l-5e 



VI - 25e + £2 



Ju 



6j;) A';. 



(39) 



Here for the case e = one finds that the massless 
state, the photon At^, no longer couples with the hid- 
den sector, while the massive mode A'^j couples with 
both the visible sector via J and with the hidden sec- 
tor via J'. We conclude, therefore, that in the ab- 
sence of the Stueckelberg mass mixing, for the case 
when only one mode is massless, there are no milli- 
charged particles coupled to the photon field. Thus 
milli-charge couplings appear in this case only when 
the Stueckelberg mixing parameter e is introduced. 
Thus for the case when only one mode is massless the 
kinetic mixing by itself docs not allow milli-charges 
but the Stueckelberg mass mixing model does. 



Appendix C: Analysis of relic density of 
milli-charged dark matter in the hidden sector 

The analysis of relic density involves the integral of 
the thermally averaged cross section from the current 
temperature to the freeze out such that [13] 



)dx. 



(40) 



Here a is the annihilation cross section for the pro- 
cess XX ~* // where / denotes a quark or a lepton, 
and V is the relative velocity of the annihilating Dirac 
fermions x ^-nd X- We use the notation, x = TjM^ 
in units where the Boltzman constant is unity, and x / 
is the value of x at the freeze out temperature. The 
thermally averaged cross section {av) is given by 



(41) 
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Using the analysis of Ref. [23| (crw) is given by 

{av) c + a) + f2ismL^B.) (42) 

where the functions /i and /2 depend on the square 
of the CM energy s and are given by 



f2{s 



is) = CfV^[8mjis + 2M^)]. 



(43) 



Here Cf = iV//(327r), Nf = (1,3) for / = 
(lepton, quark), and are as defined in Ref. [20l | 
and contain the photon, the Z and the Z' poles, the 
latter augmented with a Breit-Wigner distribution for 



thermal averaging, as the largest contribution to (av) 
for 4M^ ^ Mf arises from integration over the Z' 
pole. In the vicinity of the Z' pole we can expand 
(crv) as follows 



2 CsM-^ + CiM-^v^ 



(44) 



(v2 - e^)2 + 
where and are given by 



en = {Ml - 4M2)/Af2 , -y^, = Tz'Mz'/M^ (45) 



and where Ci{i = 1 — 4) can be read off from the 
expansion of s in Eq. and Eq. Using the 

technique of Ref. the integral J(x/) is given by 



where 



J{xf) = XfCi+ 3x^02 + 



—h + Mldih+h) 



/2 - — 



1 f 



4 + 7I 



7i?. Jo 



l/2g-y 


arctan 


4x/y - efl. 


+ arctan 



















(46) 

(47) 
(48) 



From the above one then obtains the relic density using the standard relations as discussed, for example, in 
Ref. H. ^ 



J 



Appendix D: The effect of 5 on the cross section 

The kinetic mixing parameter 5 enters Eq. (|3ip via the hidden sector where 

= x7M4'^; + 4^M + ^?^ax- 

In the limit e, (5 ^ 1, the couplings to leading order are 



(49) 



e|' = gxQx 
e| = gxQx 



1 



gxQx 



7^l2 
7^l3 



s 

6 

7T^ 



7^2l 
7^22 
7^23 



gxQx, 

egxQx sin 9 
-e-gxQx cos 6* 



1 



(50) 
(51) 
(52) 



For the purpose of illustration in this Appendix, and only here, we will work in the limit m/ ^ , w ~ and 
r^' so that 



27r 



(53) 



12 



where for 4M^ > M|, 



- "^.^ + 4M2 _ Af|/ (54) 



Defining x = (4Af ^ — Af|, ) , crw can be further simpUfied as 

2b c 
av ~ a H 1 — TT. 

where a, 6, c are functions of ?/ = (1 + Sjt) and are given by 

X J / 2 , 2 \ 



(55) 



2tt 
MlNf 



2tt 

2 ) 

(56) 



27r 

The 77 dependence can then be read off so that 0(77) — rfa{\),h{rf) = rjb{l),c{r]) = c(l). The cross section 
described by Eq. ([55]) has a dip when Xmin = —c/b which gives {(Jv)min = a — b^ /c. Now Xmin and {(Jv)min 
have an rj dependence of the form Xmm(ji) — Xmin{^) /rj, and (crw)mm(?7) = if'{'yv)min{^) which imphes that for 
the case 77 > 1 a larger rj shrinks Xmin and moves the dip closer to the resonance point. At the same time it 
also increases {(Jv)min- More interestingly, a negative t] will switch the dip to the other side of the resonance. 
These effects are illustrated in the analysis of Fig. ([6]). The analysis above explains why the inclusion of kinetic 
mixing enhances the region where the relic density constraint is satisfied. 
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